Screening Modifications of Gravity through Disformally Coupled Fields 
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It is shown that extensions to General Relativity, which introduce a strongly coupled scalar field, 
can be viable if the interaction has a non-conformal form. Such disformal coupling depends upon 
the gradients of the scalar field. Thus, if the field is locally static and smooth, the coupling becomes 
invisible in the Solar System: this is the disformal screening mechanism. A cosmological model is 
considered where the disformal coupling triggers the onset of accelerated expansion after a scaling 
matter era, giving a good fit to a wide range of background observational data. Moreover, the 
interaction leaves signatures in the formation of large-scale structure that can be used to probe such 
couplings. 
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In the standard ACDM model of cosmology, the uni- 
verse at the present day appears to be extremely fine 
tuned. The energy scale of the A component is about 
10~ 30 times the most naive expectations of theory, and 
this A has begun dominating the universal energy budget 
at a redshift which, compared to the redshift to the for- 
mation of simplest elements, is a fraction of about 10 -11 . 
In attempts to avoid these fine tuning problems, A is of- 
ten generalized to a dynamical scalar field, whose time 
evolution could more naturally result in the observed en- 
ergy density today [lj]. 

High energy physics generically predicts an interaction 
between scalar degree of freedom and other forms of mat- 
ter, which in turn could help to explain why the field 
becomes dynamically important at the present epoch. 
There are myriad variations of such models, but in all 
of them the coupling can be effectively described by a 
field-dependent mass of the dark matter particle. Those 
Yukawa-type couplings can be motivated by a conformal 
relation to scalar-tensor theories, which includes also the 
f(R) class of modified gravity [2j. 

However, for any other type of gravity modification, 
the relation between the matter and gravitational metric 
will be non-conformal. This can also be motivated e.g. 
in a DBI type scenario where matter is allowed to enter 
the additional dimensions |3j- When given by a scalar 
field </>, the disformal relation can be parametrized as 



9^u = C(<p)g^ + D 



(1) 



where commas denote partial derivatives. Considering 
the most general physical case, Bekenstein Q argued that 
both functions C and D may also depend upon (d(f>) 2 , 
but we will focus on the simpler case here. Previous ap- 
plications of such a relation to cosmology include vary- 
speed of light theories [5j, inflation [6j, dark energy 
§|, gravitational alternatives to 0, and extensions 
11] dark matter. The disformal generalization of cou- 
pled quintessence here introduced is a simple set-up that 
is useful to study generic features of the relation ([T]) in 
different scenarios. 



The disformal matter coupling. Our aim is to explore 
the novel features from the disformal coupling in the min- 
imal setting where gravity is Einstein's and the scalar 
field is coupled to a single matter species. The two met- 
rics enter into the action for gravity and the coupled 
scalar-matter system in mutually exclusive sectors 
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where the matter Lagrangian is constructed using g^ v 
from Eq. ([T]). The stress energy tensors definition 
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ensures that Einstein field equations have the usual form 
= SttGT^ . However, the covariant conservation of 
energy momentum does not hold for the coupled compo- 
nents separately. Instead, we obtain that 
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The coupling will then generically involve second deriva- 
tives, which entail the distortion of causal structure. 

For a point particle, and taking into account the cor- 
rect weight of the delta function, we have 
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From the point of view of the physical frame, the proper 
time the particle experiences is dilated by the conformal 
factor C. In addition, the disformal factor D gives a 
direction-dependent effect proportional to the projection 
of the four-velocity along the gradient of the field: 
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Extremising the proper time of the particle along its path 
shows that it follows the disformal geodesies: 



r^x? = 
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where T^p has a rather complicated form involving sec- 
ond derivatives of the field [12J . 

It is possible to derive the field equation without mat- 
ter energy-momentum derivatives using projections of 
Eq. © EJ: 

MTVySA + c ° 2DX QvJTg + V = , (7) 



where M"" = C^gT + 
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and V = £<p,<p + 2XC,^.x$- A canonical field is assumed 
in the following, = X — V . 

Equation ([7]) is a quasi-linear, diagonal, second or- 
der, partial differential equation. Its hyperbolic char- 
acter depends on the signature of the tensor Ai^, which 
involves the coupled matter energy-momentum tensor. 
For a perfect fluid, in coordinates comoving with it, 

^ v = 5 u ~ c-2DX dia g(~ P,P,P,P)- Positive energy 
density keeps the correct sign of the time derivative term 
if D > 0. However, a large pressure can flip the sign of the 
spatial derivatives coefficient, introducing an instability. 
The present analysis focuses on non-relativistic environ- 
ments, and hence q^dx ^ ^ w ^ ^ e further assumed. 
Future work will address the effects of pressure, includ- 
ing the circumstances under which the stability condition 
can break down dynamically fl2j ]. 

An example cosmological model. Let us consider an 
application where the field acts as quintessence and the 
disformal coupling is used to trigger cosmic acceleration. 
The Friedmann equations have the usual form 
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but the conservation equations for matter and the scalar 
field have to be computed from ©, ©: 



p + 3Hp = Qo4>, 



+ 3H<j> + V' = -Qo , (8) 



were the background order coupling factor reads 
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after solving away the higher derivatives. In the following 
we restrict to flat space, K = 0. 

To study the dynamics, we specify an exponential 
parametrization for the disformal relation and the scalar 
field potential 

C = C e^ /M p , D = D oe W-*°)/- M * , V = T/ o e~ 70/JUp 



with M p = (8-7rG) -1 / 2 . Besides being motivated from 
some high energy scenarios, the exponential forms fa- 
cilitate the choice of natural scales for the constant 
prefactors by shifting the zero point of the field (e.g. 
D ~ Mp 4 , Vq ~ M^, Co dimensionless) . Furthermore, 
these forms allow a convenient exploration of the phase 
space of the system. In addition to the previously stud- 
ied fixed points [3, [T3| , we find only one new, a disformal 
scaling solution that is not an attractor [12 1. 



Here we present numerical results for an example 
model where the relation (JT]) is purely disformal (Co = 1, 
a = 0) and only affects the dark matter component. We 
include radiation and baryons to consider the full realis- 
tic universe model. In this scenario, the early evolution 
is as in the usual exponential quintessence model, where 
the slope of the potential 7 has a lower limit due to the 
presence of early dark energy f2 d e = 3(1 + w)/j 2 
The new features appear when the disformal factor Dcp 2 
grows towards order one. Then the clocks that tick for 
dark matter, goo = — 1 + Dcj) 2 , slow down and make the 
effective equation of state for dark matter approach mi- 
nus unity asymptotically. The field also begins to freeze 
to avoid a singularity in the effective metric <? MI ,, and the 
universe enters into a de Sitter stage. This natural re- 
sistance to pathology was also observed in the disformal 
self-coupling scenario 0, IH ■ 

Thus, the disformal coupling provides a mechanism 
that triggers the transition to an accelerated expansion. 
The relatively steeper the slope of the disformal function 
is, i.e. the higher the ratio /3/7, the faster the transition 
happens, as seen in FIG[T] This transition also produces 
a short "bump" in the equation of state, which may have 
interesting observational consequences. We performed a 
full background MCMC analysis with a modified version 
of CMBEasy [3| using the Union2 Supernovae compi- 
lation [lj|, WiggleZ baryon acoustic scale data [l6|, cos- 
mic microwave background angular scale [17[ and bounds 
on early dark energy [lij]. The obtained constraints are 
shown in FIG [2] We see that for steep slopes 7 and 
/3, the background evolution becomes increasingly sim- 
ilar to ACDM. At this level there are no higher bounds 
on these parameters, and the model is completely viable 
with Xdist = 538.79 versus Xacdm = 538.91 (best fit 
WMAP7 parameters). However, the model is essentially 
different from ACDM, as is quite obvious when one looks 
at the effective dark matter equation of state in FIGHJ 

Cosmological Perturbations. A more realistic descrip- 
tion requires considering cosmological perturbations. In 
the Newtonian gauge, the linearized field equation is 



3ff <ty + -r + V 



SQ - 2$ (Qo + V) 
-3*), (10) 



FIG. 1: Equation of state for the field (red) and coupled mat- 
ter (blue) for different choices of the coupling slope ft. High 
values of /J/7 (solid, dashed) give a good fit to observations, 
while low values (dotted) do not produce enough acceleration. 



while the perturbed continuity and Euler equations for 
coupled dark matter are 
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The general coupling perturbation 5Q is a much more 
cumbersome combination of the fluid and field pertur- 
bations than in the purely conformal coupling 5Q( C > — 
ilog(C)'p5 + ilog(C)>^0- 

To extract the most relevant new features by analytic 
means, we shall consider the subhorizon approximation. 
In the small scale limit, taking into account only the mat- 
ter perturbations and the gradients of the field, there is 
a simple expression for the perturbed interaction SQ. In 
this Newtonian limit, we further relate the field gradient 
to the matter perturbation through the field equation 
([TO)) , which yields the simple expression SQ^ = Qo$- 
Combining the dark matter equations (fTTj) and (IT2l to- 
gether with the usual Poisson equation, we obtain the 
evolution of the coupled dark matter overdensity 
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The Hubble friction is accompanied with an additional 
term due to the evolution of the field, and the source term 
is modulated. The last effect is captured by defining an 
effective gravitational constant G e f / that determines the 
clustering of dark matter particles on subhorizon scales 

*ff _ 1 = . (14) 



G 



G ' AirGp 2 ' 

This approximation has the same expression as the sim- 
ple conformal case, although with a significantly different 
functional form of the coupling ©. 



FIG. 2: Marginalized one and two-sigma regions obtained 
from Supernovae (Blue), BAO (Green), CMB angular scale 
+ early dark energy bounds (Orange), and combined con- 
straints. All contours included a prior on Ho from the HST 
[2l| and QbHg from Big Bang Nucleosynthesis [22T ], 



For our example model, the late time dependence dur- 
ing dark energy domination produces a large enhance- 
ment of the matter growth, SG e ff/G ~ ("fV/p) 2 3> 1, 
as 7 > 10 is required to avoid the effects of early dark 
energy. Such behavior is in tension with large scale struc- 
ture observations, and also occurs in conformally coupled 
models that attempt to address the coincidence problem 
[l9j |. This discrepancy can be alleviated with a differ- 
ent choice of the functions C,D,V. One possibility is to 
introduce a modulation D{4>) — > f(<j>)D(<fr) to make Qo 
small enough after the field enters the slow roll phase. 
This modification can render SG e ff arbitrarily small, ex- 
cept for a relatively short time around the transition (see 
FIGU]). A different solution to control perturbations can 
be achieved by allowing the field itself to live in a disfor- 
mal metric, as in the uncoupled models studied in Refs. 
0, H| ■ Viable variations might be exploited to alleviate 
the claimed problems of ACDM with small scale struc- 
ture formation [20j |. 

The disformal screening mechanism. Finally, we con- 
sider the possibility of extending the disformal coupling 
to visible matter. Due to the stringent bounds on equiv- 
alence principle violations 



23j, some sort of screening 



mechanism is necessary to hide the coupling in dense 
environments such as the Solar System. The effects of 
disformal couplings to the chameleon screening were re- 
cently investigated by Noller 24j, who correctly noted 
that the disformal contribution to the conservation equa- 
tions vanishes for static, pressureless configurations. This 
is obvious from ((5J, since only the T 00 component is 
nonzero for dust, and when contracting the field deriva- 
tives with the stress tensor a non-vanishing result re- 
quires time evolution of the scalar field. 
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Addressing the effects of disformal couplings therefore 
requires studying the field dynamics in high density, non- 
rclativistic environments. This regime can be explored 
using the general scalar field equation for a station- 
ary density distribution p(x) in the limit p — > oo, and 
neglecting the remaining spacetime curvature 31]. The 
same result follows from taking the limit p ^> C/D, cb 2 in 



2£T \ C 2D 



2M„ 



(15) 



where the first equality is general and the second applies 
to our example model. The above expression departs sub- 
stantially from the simple conformal coupling, for which 
the p — > oo limit is ill-defined. Spatial derivatives become 
irrelevant, as they are suppressed by a p/ p factor w.r.t. 
time derivatives. More importantly, the above equation 
becomes independent of the local energy density, making 
the field evolution insensitive to the presence of massive 
bodies. As the field rolls homogeneously, spatial gradi- 
ents between separate objects, which would give rise to 
the scalar force, do not form. 

In the purely disformal case with exponential D, equa- 
tion (fT5|) can be integrated directly 
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In this solution, the field time variation is approximately 
constant while 4 <C 4^ cx </>(0) _1 and slows down after- 
wards as cx 1/4. Since the coupling is proportional to <f>, 
stronger couplings decay earlier. Furthermore, as V is a 
decreasing function of the field, the interaction density 
p = — {(f) + V')<fi has at most cosmological values and its 
effect on dense systems is highly suppressed. Assuming 
(ft 2 V ~ po, Dpo > 1, where po is the average cosmic 
density, typical mass increase rates M/M are as small as 
~ 10~ 6 /Gy for the interstellar medium and ~ 10~ 29 /Gy 
for the average Earth density. The induced mass vari- 
ation would have consequences e.g. in the Earth-Moon 
system, altering the orbits in a way degenerate with a 
time evolution of Newton's constant. However, the pre- 
dicted rates are way below current lunar laser ranging 
bounds G/G < 10" 3 /Gy HH. 

New, potential signatures may be found in the presence 
of matter velocity flows, non negligible pressure, field 
gradients of cosmological origin, or in the mildly non- 
linear regime, when the screening starts taking place [lij ]. 
Other proposed tests involve high-precision, low-energy 
photon experiments j2(| and effects on the chemical po- 
tential of the baryon-photon fluid [27| . 

This provides a novel disformal screening mechanism, 
which is distinct from the Vainhstein [28| , chameleon [29[ 
and symmetron [30j | mechanisms, which respectively rely 
on the nonlinearity of the field derivatives and the de- 
pendence of the field mass and coupling with the ambi- 
ent density. Our mechanism relies on the existence of a 



well defined p — > oo, non-relativistic limit in the scalar 
field equation, given by equation (|15[) . for which the field 
evolution is independent of the matter energy density. If 
the conformal part C is negligible, only a friction term 
remains and the field coupling density ^ is a decreas- 
ing function of time. As it evolves below its cosmological 
value (provided V < and D'/D > 0), the effects of 
the coupling are suppressed by a factor ~ Po/p and the 
theory is consistent with precision gravity tests. 

In summary, the idea of disformal coupling survives 
scrutiny with so high flying colours that the reader must 
be suspicious too. The theory has a rich structure and 
many interesting applications. At cosmological scales, it 
can mimic closely the ACDM evolution, providing means 
to address the coincidence and cosmological constant 
problems. At cluster and galactic scales, new effects are 
to be expected due to the novel manifestation of the fifth 
force that could help to cure the discrepancies that the 
standard model has with large scale structure observa- 
tions. Finally the disformal screening mechanism oper- 
ating on very small scales opens a completely new avenue 
for fundamental scalar fields strongly coupled to matter 
into our reality. 

We thank L. Amendola and J. Garcfa-Bellido for dis- 
cussions. TK and DFM are supported by the Norwegian 
research council. MZ is supported by MICINN through 
BES-2008-009090. 



[1 

[2 

[3 

[4 
[5 
[6 

[7: 

[8 

[9 
[10 

[11 

[12 

[13 
[14 

[is; 

[16 

[17 
[is; 

[19 



E. J. Copeland, M. Sami, and S. Tsujikawa, 
Int.J.Mod.Phys. D15, 1753 (2006), hep-th/0603057. 
T. Clifton, P. G. Ferreira, A. Padilla, and C. Skordis, 
Phys.Rept. 513, 1 (2012). 

C. de Rham and A. J. Tolley, JCAP 1005, 015 (2010), 
1003.5917. 

J. D. Bekenstein, Phys.Rev. D48, 3641 (1993). 

J. Magueijo, Rept.Prog.Phys. 66, 2025 (2003). 

N. Kaloper, Phys.Lett. B583, 1 (2004), hep-ph/0312002. 

T. S. Koivisto (2008), 0811.1957. 

M. Zumalacarregui, T. Koivisto, D. Mota, and P. Ruiz- 

Lapuente, JCAP 1005, 038 (2010), 1004.2684. 

J. D. Bekenstein, Phys.Rev. D70, 083509 (2004). 

M. Milgrom, Phys.Rev. D80, 123536 (2009), 0912.0790. 

D. Bettoni, S. Liberati, and L. Sindoni, JCAP 1111, 007 

(2011) , 1108.1728. 

M. Zumalacarregui, T. S. Koivisto, and D. F. Mota 

(2012) , see also suplemental material, 1210.8016. 
L. Amendola, Phys.Rev. D62, 043511 (2000). 
M. Doran, JCAP 0510, 011 (2005). 

R. Amanullah, C. Lidman, D. Rubin, G. Aldering, 
P. Astier, et al., Astrophys.J. 716, 712 (2010), 1004.1711. 

C. Blake, E. Kazin, F. Beutler, T. Davis, D. Parkinson, 
et al., Mon.Not.Roy.Astron.Soc. 418, 1707 (2011). 

D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, 
M. Nolta, et al., Astrophys.J.Suppl. 192, 16 (2011). 

C. L. Reichardt, R. de Putter, O. Zahn, and Z. Hou, 

Astrophys.J. 749, L9 (2012), 1110.5328. 

T. Koivisto, Phys.Rev. D72, 043516 (2005). 



■5 



[20] L. Perivolaropoulos (2011), 1104.0539. 

[21] A. G. Riess, L. Macri, S. Casertano, H. Lampeitl, H. C. 

Ferguson, et al., Astrophys.J. 730, 119 (2011). 
[22] K. Nakamura et al. (Particle Data Group), J.Phys.G 

G37, 075021 (2010). 
[23] C. M. Will, Living Rev.Rel. 9, 3 (2005). 
[24] J. Noller (2012), 1203.6639. 

[25] J. G. Williams, S. G. Turyshev, and D. H. Boggs, 

Phys.Rev.Lett. 93, 261101 (2004), gr-qc/0411113. 
[26] P. Brax, C. Burrage, and A.-C. Davis (2012), 1206.1809. 



[27] C. van de Brack and G. Sculthorpe (2012), 1210.2168. 
[28] A. Vainshtein, Phys.Lett. B39, 393 (1972). 
[29] J. Khoury and A. Weltman, Phys.Rev.Lett. 93, 171104 
(2004). 

[30] K. Hinterbichler and J. Khoury, Phys.Rev.Lett. 104, 

231301 (2010), 1001.4525. 
[31] More precisely, neglecting p/p, (C — 2DX)/Dp, Poo^n, 

(C - 2DX)V/(Dp), ptfJip^ 2 ) and assuming Dp < 

C - 2DX. 



